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Sine-Gordon  Phase  Space  -  Geometry  and  Instabilities 

by 


Nick  Ercolani 


Abstract: 

The  Sine-Gordon  Hamiltonian  system 


on  the  phase  space  of  functions  u  =  [v(xj]  e  CL  x  CL'  Peri°d^c  i-n  x  of 
period  L,  is  completely  integrable.  The  complete  involutive  set  of 
integrals  is  encoded  in  the  coefficents  of  an  entire  function  of  X  e  C 

A  (u;  X) 

which  depends  on  u.  (A  is  the  Floquet  discriminant  of  a  X-eigenvalue 
problem  whose  coefficients  involve  u.) 

The  modulational  instabilities  of  Sine-Gordon  are  associated  to 
the  level  sets  of 

A:  c*  x  c*  - ♦  {entire  functions} 

L  Li 

which  lie  over  critical  values  of  A  (={A(X)|A±2  has  complex  double 
roots) . 

Such  a  level  set  is  stratified  by  a  finite  number  of  invariant 
manifolds.  We  show  that  the  topology  of  these  manifolds,  how  they  are 
connected  to  one  another,  ar.d  how  this  topology  governs  the  Lyapunov 
exponents  of  associated  ir.st  abilities  can  all  be  determined  from  the 
function  A (X) . 
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Spatial  Coherence  and  Temporal  Chaos  in  Near-Integrable  PDE's 


by 

David  W.  McLaughlin 

Abstract: 

The  damped  and  driven  pendula  chain, 

utt  -  u  +  sin  u  =  €  [-a  u t  +  T  sin  ut] , 
is  studied  as  a  prototype  model  of  coherence  and  chaos  in  a  near 
integrable  PDE.  Both  numerical  and  theoretical  studies  are  summarized. 
The  two  main  points  are  that  the  nearby  integrable  system  provides  (i) 
homoclinic  orbits  which  act  as  sources  of  sensitivity  (chaos)  and  (ii) 
candidates  for  approximate  coordinates  of  the  attractor.  In  addition, 
techniques  founded  on  the  nearby  integrable  system  provide  direct 
numerical  checks  on  the  validity  and  accuracy  of  these  two  points. 
These  numerical  and  theoretical  techniques  are  used  to  study  in  detail 
a  route  to  chaos  which  involves  the  following  characteristics: 

(i)  temporal  -  one  frequency  -*  two  frequencies  -*  chaos; 

(ii)  spatial  -  zero  -»  one  ->  two  localized  excitations; 

( i i i )  symmetry  changes  and  pattern  competition; 

(iv)  low  dimensional,  chaotic  attractors; 

(v)  temporal  intermittancy ; 

(vi)  homoclinic  crossings; 

(vii)  interactions  and  transitions  between  localized  and  extended 
states 


References: 


1)  N.  Ercolani,  G.  Forest,  D.  McLaughlin  (i)  Lect.  Appl.  Math.  23, 
149-165  (1986);  (ii)  Physica  18D,  472-474  (1986);  (iii)  "Homoclinic 
Orbits  for  the  Periodic  Sine-Gordon  Equation",  submitted,  Physica  D; 

(iv)  "Geometry  of  the  Modulational  Instability"  Parts  I  and  II 
(Preprints,  University  of  Arizona). 

2)  A.  Bishop,  M.  G.  Eorest,  D.W.  McLaughlin,  E.  Overman,  Physica 
D.  1 1,  293-328  (1986) . 

3)  (i)  A.  Bishop,  D.W.  McLaughlin,  E.  Overman,  to  appear,  Proc.  of 
Conf .  on  Solitons,  Tech.  Univ.  Denmark,  Ed.  by  P.  Christensen  (1987) ; 
(ii)  Physica  19D,  1-41  (1986) . 

4)  A.  Bishop,  et.al.  Physica  TD,  759-779  (1983). 

5)  N.  Ercolani  and  M.  G.  Forest,  Comm.  Math.  Phys.  99,  1-45 
(1985)  . 

6)  H.  McKean,  Comm.  Pure.  Appl.  Math  3j£,  (1981). 

7)  Physica  7D  -  The  entire  volume;  (ii)  Physica  23D  -  The  entire 
volume. 
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Correlations  Between  the  Perturbed  Sine-Gordon  Equation 
and  Finite  Modal  Equations 

by 

Greg  Forest 

Abstract: 

In  this  lecture  we  describe:  i)  numerical  results  on  the 
bifurcations  of  the  damped,  periodically  forced  Sine-Gordon  equation 
with  periodic  boundary  conditions,  in  a  finely  tuned  parameter  range; 
ii)  an  interpretation  of  the  spatial  and  temperal  bifurcation 
structures  of  this  perturbed  integrable  system  with  regard  to  the  exact 
structure  of  the  Sine-Gordon  phase  space;  iii)  a  model  dynamical 
systems  problem,  which  is  itself  a  perturbed  integrable  Hamiltonian 
system,  derived  from  the  perturbed  Sine-Gordon  equation  by  a  finite 
mode  truncation  in  the  nonlinear  Schrodinger  limit;  and  iv)  the 
bifurcations  to  chaos  in  the  four  dimensional  model  problem. 

In  particular,  we  focus  on  a  likely  source  of  chaos  in  both  the 
o.d.e.  and  p.d.e.  systems:  the  existence  of  homoclinic  orbits  in  the 
unperturbed  integrable  phase  space  and  the  continuation  of  these 
homoclinic  structures  in  the  perturbed  problem.  Finally  we  numerically 
correlate  the  homoclinic  crossings  in  the  chaotic  dynamics  of  the  full 
and  reduced  problems. 

These  experimental  results  provide  physical  intuition  about  the 
coexistence  of  simple  coherent  spatial  structures  and  temporal  chacs, 
and  set  the  stage  for  carrying  out  the  rigorous  mathematical  analysis 
to  support  the  numerical  work. 


My  collaborators  are  Nick  Ercolani  and  Dave  McLaughlin  on  the 
theoretical  aspects,  Alan  Bishop  on  the  formulation  of  the  experiments 
and  Ed  Overman,  Yannis  Kevrekidis,  Mike  Jolly,  Mac  Hyman,  and  Randy 
Flesch  on  the  numerical  studies. 

References : 

1.  "A  quasi-periodic  route  to  chaos  in  a  near-integrable  p.d.e.,"  A. 
Bishop,  G. Forest,  D.  McLaughlin,  and  E.  Overman,  Physica  23D 
(1986),  293-328,  and  references  therein. 

2.  "Correlations  between  the  perturbed  Sine-Gordon  equations  and 
finite  model  equations,"  A.  Bishop,  R.  Flesch,  G.  Forest,  D. 
McLaughlin,  and  E.  Overman,  preprint,  May,  1987. 

3.  "Geometry  of  the  modulational  instability:  local  results"; 
"Geometry  of  the  modulational  instability:  gobal  results"? 
"Homoclinic  orbits  in  the  periodic  Sine-Gordon  equation",  by  N. 
Ercolani,  G.  Forest,  and  D.  McLaughlin,  preprints,  1987. 
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An  Infinite  Dimensional  KAM  -  Theorem 

by 


J.  Poschel 

Abstract; 

We  consider  a  d-dimensional  lattice  A  =  Zd  of  harmonic 

oscillators,  whose  frequencies  o^,  i  e  A,  are  considered  as  parameters 

which  may  be  adjusted  if  necessary.  The  unperturbed  situation  is 

described  by  the  Hamiltorian 

N  =  (u ,  y )  =  .2.  u.y. 

v  111  leA  in 

in  the  phase  space  I7A  x  RA .  We  consider  a  perturbation  H  =  N  +  aP, 

where  P  has  a  spatial  structure  (analogous  to  a  Fourier  series 

expansion) :  P  =  2  Pa  , 

Aed  A 

where  s<  is  a  system  of  finite  subsets  of  A,  and  PA  "lives  on  A".  The 
size  of  the  perturbation  is  expressed  in  terms  of  a  measure  u  on  A ,  fcr 
example; 

(*)  M (A)  =  I  | i  | 

i€A 

We  then  require  (roughly  speaking) 

|PAI  -e-a"<A>, 

with  some  suitable  a  >  o. 

We  are  able  to  shew  that  then,  for  sufficiently  small  a,  the 
perturbed  system  possesses  a  Cantor  set  of  infinite  dimensional, 
invariant  smooth  tori  with  linear  flow  on  them.  The  frequencies  u  cf 
these  tori  satisfy  the  snail  Jivision  condition 
(**)  |  (k,  «)  I  >  -r  6(  !>:!  ,  '  •  ••.:>)  ,  0  Jt  k  €  ZA  , 


6 


where 


<k>  =  y (supp  k) 

...  1+6.. 

5 (t)  =  e"C/iog  ,  t  l  tQ/  a  >  0. 

We  remark  that  for  example  with  (*),  condition  (**)  holds  on  a  set  of 
large  (Gaussian)  measure  in  u-space  R^ . 

«• 

The  above  generalises  and  improves  a  result  by  Frohlich,  Spencer 
and  Wayne  (/.  Stat.  Phys  .  1986.42) . 


Principle  of  Spatial  Averaging  and  Inertial  Manifolds 


by 


Ihe  key  idea  in  the  lecture  is  to  introduce  the  Principle  of 

Spatial  Averaging  and  to  show  this  principle  can  be  used  to  show  the 

existence  of  inertial  manifolds.  This  theory  does  apply  in  the  case 

3  3 

where  (1)  is  dissipative,  f  €  c  ,  and  n  =  [0,2tt]  ,  where  the  boundary 
condition  are  of  Dirichlet,  Neumann  or  periodic,  types. 
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Integral  Manifolds:  Inertial  Manifolds 

by 

Ciprian  Foias 

Abstract: 

Many  partial  differential  equations  describing  dissipative 
phenomena  have,  when  supplemented  with  appropriate  boundary  conditions, 
the  following  form 

(1)  ut  +  Au  +  R(u)  =  0 

where  A  is  a  positive  operator  on  a  Hilbert  space  with  a  compact 
inverse  A  1  and  R(u)  is  a  differentiable  operator  from  the  domain  2(A) 
of  A  into  3f,  satisfying 

(2)  |R  (u)v|  i  c | Au  j  | A^v |  (u , v€2 (A) ) 

where  c  >  0,  0  i  p  <  1/2  are  constants  and  |* |  denotes  the  norm  in  % . 
An  inertial  manifold  of  (1)  is  a  finite  dimensional  (Lipschitz) 
manifold  J  in  2(A),  invariant  to  (1)  (i.e.  the  solutions  of  (1)  which 

start  in  J  remain  in  M )  and  exponentially  attracting  all  bounded  sets 
in  3f  [1).  A  natural  way  to  construct  such  a  manifold  is  to  choose  a 
large  sphere  r  in  the  linear  space  generated  by  the  eigenvectors  of  A 
corresponding  to  the  first  m  distinct  eigenvalues  <  A <♦••<  Am  of 
A,  and  to  construct  the  integral  manifold  2  of  (1)  determined  by  T;  of 
course  T  must  not  be  characteristic,  i.e.  Au  +  R(u)  should  not  be 
tangent  to  T  at  any  u  €  r.  Finally,  one  shows  that  the  closure  i  of  - 
in  2(A)  is  an  inertial  manifold.  This  method  was  developed  in  [2]  fcr 
several  remarkable  partial  differential  equations.  In  [2]  a  central 
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role  is  played  by  the  spectral  blocking  property  of  2.  Namely,  for  u  e 

2  let  A (u)  denote  the  maximum  of  (Ag,g)  over  all  vectors  g  of  length  < 

1  and  tangent  to  2  at  u.  Assume  that  the  spectral  gap  between  and 

the  next  eigenvalue  A^+1  of  A  is  large  enough;  more  precisely  that 

,,,  Am+l"Am  .  _,Am+l+Anu0  , .  , 

(3)  - - -  >  C( - - - )K  sup  |au|. 

u€2 

Then  if  for  all  u  e  r  we  have 


A  (u)  i 


A  .  +A 
m+1  m 


this  last  inequality  holds  for  all  u  €  I.  a  sketch  of  the  proof  of 
this  property  was  given  in  the  lecture,  some  remarks  and  conjectures 
in  connection  with  the  above  topics  were  also  presented;  particularly  a 
candidate  for  an  inertial  manifold  of  the  2D  Navier-Stokes  equations 
was  introduced. 

[1]  C.  Foias,  G.R.  Sell,  R.  Temam:  Inertial  manifolds  for  nonlinear 
evolutionary  equations,  IMA  preprint  series  #234  (U.  of  Minnesota, 

March  1986) . 

[2]  P.  Constantin,  C.  Foias,  B.  Nicolaenko,  R.  Temam;  Nouveaux 

1  t  , 

resultats  sur  les  varietes  inertielles  pour  les  equations 

» 

differentielles  dissipatives,  C.R.  Acad.  Sci.  Paris,  302,  Serie  I, 

1986,  375-378;  Integral  manifolds  and  inertial  manifolds  for 
dissipative  partial  differential  equations  (in  preparation). 
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Inertial  Manifolds  For  Dissipative  Perturbed  Hamiltonnian  Systems 


Basil  Nicolaenko 


We  present  an  outline  of  a  general  theory  of  Inertial  Manifolds 
for  damped  Hyperbolic  systems  and  dissipatively  perturbed  Hamiltonian 
Systems.  A  prior  lack  of  compactness  of  the  relevant  semi-groups 
raises  the  usual  problems. 

We  introduce  appropriately  modified  Hamiltonians  to  penorm  the 
problem.  The  "strong  squeezing"  and  "cone  invariance"  properties  are 
still  essential  in  constructing  Inertial  Manifolds.  In  contrast  with 
previous  work  on  parabolic  dissipative  systems,  these  "cone  propertie 
now  only  make  sense  with  respect  to  the  metric  (equivalant  norm) 
induced  by  the  Hamiltonian. 

We  give  an  explicit  construction  for  a  system  of  perturbed 
conservation  laws  modeling  chaotic  vapor-liquid  phase  changes  in  a 
compressible  flow  with  non-convex  Van-der-Waals  equations  of  state. 
Computer  movies  show  complex  temporal  interactions  between  large  seal 
spatial  structures;  this  suggests  a  reduced  set  of  local  coordinates 
the  corresponding  inertial  manifolds. 


1.  "Integral  and  Inertial  Manifolds",  P.  Constantin,  C.  Foias,  B. 

Nicolaenko,  R.  Terr,  submitted  to  Springer  Verlag  Lecture  Notes 
in  Mathematics. 
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2.  "Inertial  Manifolds  for  Dissipative  Perturbed  Hamiltonian 

Systems",  B.  Nicolaenko,  (future  M.S.I.  Report);  also  submitted  to 
"A.M.S.  Contemporary  Mathematics"  Series. 


Differentiable  Foliation  in  Infinite  Dimensional  Systems 


Shui-Nee  Chow 


We  consider  abstract  evolution  equations  in  a  Banch  space  Z. 

x  =  A  x  +  F (x)  ,  x  e  z 

O 

where  A  generates  a  C  semigroup  of  linear  operators  on  Z  and  F  is  a 
nonlinear  function.  Under  some  general  conditions,  we  show  that  if  F 

3c  CL 

is  C  ,  then  there  exist  a  C  ,  o  <  a  <  1,  foliation  near  x  =  0. 
Furthermore,  each  fiber  is  C  .  The  relation  between  foliations  and 
inertial  manifolds  is  discussed.  Generalization  to  compact  invariant 
sets  is  also  considered. 


[1]  S.N.  Chow,  K.  Lu,  and  Xiao  Biao  Lin,  Differentiable  Foliations  in 
Infinite  Dimensional  Systems.  In  preparation. 


Remarks  on  the  Ergodic  Theory  of  Foliations  and  Attractors 


by 


M.  Brin 


The  main  purpose  of  this  talk  is  to  present  some  general 
principles  of  studying  finite  dimensional  dynamical  systems  with 


attractors  by  means  of  ergodic  theory  and  to  indicate  a  possibility  of 
generalzing  these  methods  for  infinite  dimensional  evolution  equations 
with  compact  attractors.  If  a  dynamical  system  has  an  attractor,  one 
usually  hopes  that  the  asymptotic  behavior  of  the  system  when  t  -»  00  is 
very  well  modeled  by  the  induced  asymptotic  behavior  on  the  attractor. 
To  make  this  work,  one  considers  the  set  of  "good"  points  on  the 
attractor,  i.e,  those  points  whose  orbits  are  correctly  distributed  on 
the  attractor  (in  the  sense  of  the  Birkhoff  ergodic  theorem)  and  for 
which  the  linearized  equation  has  correct  asymptotic  properties  (in  the 
sense  of  the  Oseledetz  multiplicative  ergodic  theorem) .  The  next  step 
is  to  construct  the  stable  manifolds  of  "good"  points  and  to  conclude 
that  any  point  in  the  stable  manifold  of  a  "good"  point  has  the  same 
asymptotic  characteristics  (for  t  -♦  +«)  as  the  "good"  point  itself. 

The  final  step  is  to  show  that  almost  every  (w.r.t.  the  Lebesgue 
measure)  point  from  the  basin  of  attraction  has  correct  asymptotic 
properties.  To  do  that  one  shows  that  the  stable  foliation  satisfies  a 
version  of  the  Fubini  theorem,  i.e.  is  absolutely  continuous. 


The  first  step  in  generalizing  this  approach  for  the  infinte 
dimensional  case  was  made  by  D.  Ruelle  who  constructed  stable  manifolds 
for  the  "good”  points  on  a  compact  attractor  (R.  Mane  proved  a  similar 
theorem  for  Banach  spaces).  Z.  Nitecki  and  myself  showed  that  in  this 
case  the  stable  foliation  is  absolutely  continuous,  which  leaves  some 
hope  that  the  above  finite  dimensional  approach  may  be  generalized  for 
some  infinite  dimensional  evolution  equations.  The  fact  that  a  Hilbert 
space  does  not  carry  any  canonical  or  natural  measure  similar  to  the 
Lebesgue  measure,  adds  another  meaning  of  the  absolute  continuity  of 
the  stable  foliation  since  it  allows  one  to  measure  the  size  of  a  set 
consisting  of  stable  leaves.  Such  a  set  is  big  if  its  cross  sections 
by  finite  dimensional  planes  transverse  to  the  foliation  have  large 
measure.  By  the  absolute  continuity,  this  property  does  not  depend  on 
the  particular  cross  section. 
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Asymptotics  of  Quasilinear  Parabolic  Equations  of  Viscoelasticity 

by 

Stuart  S .  Antman 

Abstract: 

The  motion  of  a  mass  point  on  a  spring  is  perhaps  the  most 
fundamental  problem  in  the  theory  of  oscillations.  If  the  force 
exerted  by  the  spring  on  the  mass  depends  nonlinearly  on  the  position 
and  velocity  of  the  mass  point,  then  the  analysis  of  the  governing 
unforced  differential  equations  is  routine.  A  spring  having  these 
properties  is  effectively  assumed  to  have  zero  mass.  If  the  spring 
itself  has  mass,  then  its  motion  is  coupled  with  that  of  the  end  mass. 
In  the  case  under  study  the  spring  is  assumed  to  be  nonlinearly 
viscoelastic,  so  that  its  motion  is  governed  by  a  third-order 
quasilinear  "parabolic-hyperbolic"  system.  The  coupled  system  is 
studied  in  the  limit  that  the  mass  density  of  the  spring  goes  to  zero. 
The  leading  term  of  the  resulting  asymptotic  expansion  of  the  solution 
satisfies  the  system  in  which  the  spring  density  is  zero.  It  is  shewn 
that  the  motion  of  the  end  mass  typically  does  not  satisfy  an  ordinary 
differential  equation  of  the  type  described  above.  The  formal 
asymptotic  expansion  of  the  solution  of  the  full  initial-boundary  value 
problem  accounting  for  an  mtial  layer,  can  be  fully  justified.  The 
justif iciation  relies  c.n  a  number  of  delicate  estimates  for  quasil ir.e  -> : 
parabolic  equations. 
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1)  S.S.  Antraan,  The  Paradoxical  Asymptotic  Status  of  Massless 
Springs,  University  of  Maryland  Technical  Report,  1987. 

2)  S.S.  Antman,  Asymptotic  Analysis  of  Quasilinear  Parabolic 
Equations  of  Viscoelasticity,  in  preparation. 
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Exponentially  Small  Splitting  of  Separatrices 


by 


J.  Marsden 

Abstract; 

Both  upper  and  lower  .estimates  are  established  for  the  separatrix 
splitting  of  rapidly  forced  systems  with  a  homoclinic  orbit.  The 
general  theory  is  applied  to  the  equation 

v  +  sin  ^  =  5  sin 

for  illustration;  in  this  example,  for  any  17,  0  <  77  <  tt/2,  and  0  <  e  < 
1/  0  i  6  i  6q/  for  6q  sufficiently  small,  the  separatrices  are  proved 
to  split  by  an  amount  is  no  more  than 

6C(„,  60)  exp  (-  ±  (l  -  ->]) 

where  C(tj,  5q)  is  a  constant  depending  on  n  and  6Q.  If  we  replace  6  by 
P  >  8,  then  we  have  the  sharper  estimate 

C2ep6e  ir'/2a  <;  splitting  distance  <  C1ep5e-ir/2t 
for  constants  C1  and  C2 .  In  particular,  in  this  latter  case,  the 
Melnikov  criterion  correctly  predicts  exponentially  small  splitting. 

The  techniques  developed  here  can  be  applied  to  estimate  the  thickness 
of  stochastic  layers  in  the  unfoldings  of  degenerate  singularities  and 
in  KAM  theory. 
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Hamiltonian  Structures  and  Invariant  Manifolds  for  Elliptic 
Equations  in  Cylindrical  Domains 

by 

Alexander  Mielke 

Abstract: 

We  consider  second  order  elliptic  systems  which  are  derived  by 
minimizing  an  energy  functional  I(u)  =  W(u,vu)dy.  For  cylindrical 
domains  ft  =  (0,T)  x  2  with  variables  (t,x),  the  corresponding 
Euler-Lagrange  equations  can  be  written  in  Hamiltonian  form,  viz. 

(*)  u  =  dH/dv  ,  v  =  -  dH/du 

where  u  =  du/dt,  v  =  dw/du,  and  H(u,v)  =  /2(u*u  -  W (u , v^u , u) }dx 
(cf. [1]) . 

On  the  other  hand  it  is  well-known  that  for  elliptic  equations  cn 
infinite  cylindrical  domains  (i.e.  f?  =  IR  x  2)  with  2  bounded,  all  small 
solutions  lie  on  a  finite  dimensional  center  manifold  CM  (cf.[2,3]). 

Now  it  can  be  shown  that  the  Hamiltonian  structure  of  (*)  carries  ever 
to  a  Hamiltonian  structure  for  the  ordinary  differential  equation 
describing  the  flow  on  CM.  However  this  structure  is  non-canonical  Joe 
to  the  curvature  of  the  CM. 

An  Application  to  Sa int-Venant 1 s  Problem  is  given  [4]. 

References: 

[1]  R.  Abraham  &  J.E.  Marsden  [1978],  "Foundations  of  Mechanics",  2nd 
edition,  Addison-Wes 1 ey ,  Reading,  MA. 

[2]  G.  Fischer  [1984],  2e- t : ursrannigf altigkeiten  bei  elliptischen 
Dif  ferentialgleichur.-:*-  n ,  Math.  Nachr.  115,  137,157. 
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[3]  A.  Mielke  [1986],  Reduction  of  quasilinear  elliptic  equations  in 
cylindrical  domains  with  applications.  Math.  Methods  in  the  Appl 


Sciences,  to  appear. 

[4]  A.  Mielke  [1987],  Saint-Venant 1 s  problem  and  semi-inverse 
solutions  in  nonlinear  elasticity.  Arch.  Rat.  Mech.  Anal., 
submitted. 
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Model  Complementarity  for  the  Two-Dimensional  Euler 
and  Navier-Stokes  Equations 

by 


Norman  J.  Zabusky 

Abstract: 

Continuum  nonlinear  dynamical  systems  in  more  than  one  space 

dimension  are  almost  all  mathematically  intractable.  For  example,  the 

generic  equations  of  fluid  dynamics,  the  incompressible  Euler  (u=C)  and 

Navier-Stokes  equations  in  two  dimensions, 

d.cj  =  +  <j  +  >£  < j  =  uAu,  A  'll  —  -o , 

t  t  y  x  x  y 

are  examples  with  no  time-dependent  analytical  solutions  if  the  initial 
condition  is  two  or  more  localized  vorticity  distributions,  <j(x,y,0)  : 

k<2>. 

A  major  goal  of  the  study  of  nonlinear  dynamical  systems  is  an 

analytical  understanding  that  will  allow  accurate  predictions  of  state 

variables  or  functions  of  state  variables  over  moderate-to-long  tires. 

To  accomplish  this  objective,  it  is  necessary  to  use  more  than  one 

model  to  represent  compl :  nentary  aspects  of  a  dynamical  process.  The- 

results  of  simulations  with  these  models  must  be  visualized  and 

quantitative  feature  extraction  must  be  performed  using  numerical 

diagnostic  algorithms.  ?h:s  permits  the  judicious  investigator  to 

"separate"  complex  inter  jet  :  :r.s  among  coherent  structures  into  simpler 

and  hopefully  analytical!/  *r actable  parts.  This  synergetic  appreaer 

1  2 

was  used  in  the  discover/  f  the  soliton.  ' 
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The  essential  problems  in  inviscid,  incompressible  two-dimensional 
hydrodynamics  were  not  posed  until  computer  simulations  showed: 
axlsymmc t r I zat t on  of  noncircular  distributions  of  vorticity?  merger  of 
like-signed  vorticity;  binding  of  opposite-signed  vorticity;  and 
entrainment  in  a  host  vortex  of  small  regions  of  irrotational  fluid  of 
opposite-signed  vorticity.  In  all  these  processes  one  also  observes 
gradient  intensif ication  of  vorticity  due  to  relative  transport.  Three 
key  questions  in  evolving  inviscid  or  nearly-inviscid  two  dimensional 
flows  are: 

(1)  What  are  the  mechanisms  by  which  smooth  vorticity 
distributions  aggregate  or  "condense"  into  near-circular 
regions  of  vorticity  (that  is:  the  mechanisms  of 
anixymmetrization  and  gradient  intensification)? 

(2)  What  are  the  mechanisms  by  which  two  like-signed  regions  of 
vorticity  merge  or  the  mechanisms  by  which  two 
opposite-signed  regions  of  vorticity  bind  or  entrain? 

(3)  How  do  the  filaments  (small  scales)  which  arise  in 

axisymmetrization,  merger  and  binding  affect  the  long  time 

evolution  of  the  large-scale  structures? 

Axisymmetrization,  merger  and  binding  may  be  considered  fundamental 

3  4 

physical  space  interactions  in  two-dimensional  turbulent  flow.  '  In 
spectral  jargon,  merger  corresponds  to  an  "upward"  energy  cascade  and 
f ilamentation  accounts  for  the  "downward"  enstrophy  cascade.  The 
f ilamentation  of  vorticity  during  axisymmetrization,  merger  and  binding 
is  best  understod  by  considering  the  local  corotating  or  cotranslating 

4 

streamfunction . 
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In  a  recent  review,  Melander,  Overman  and  Zabusky  discussed 

computational  and  mathematical  vortex  dynamics,  emphasizing 

two-dimensional  aspects.  Here  we  focus  briefly  on  complementary  models 

and  numerical  diagnostics  in  this  field.  Beginning  in  the  early  1970's 

finite-difference,  spectral,  pseudospectral  and  vortex-in-cell  models 

and  algorithms  became  very  popular.  The  first  three  must  include  a 

dissipative  mechanism  to  resolve  small-scale  structures.  However,  if 

gradient-scale  lengths  form  which  are  smaller  than  five  computational 

zones  then  significant  truncation  and  aliasing  errors  arise.  Thus, 

none  of  these  schemes  are  capable  of  a  general  study  of  the  Euler 

equations  or  even  of  "very  high"  Reynolds  number  flows  at  "moderate" 

times,  despite  the  increasing  availability  of  supercomputer  resources. 

However,  aspects  of  these  flows  can  be  better  studied  with  the  contour 

dynamical  algorithms6'7'8,9'10  and  the  moment  model.11  Contour 

12  . 

dynamics,  a  generalization  of  the  "vaterbag"  model  is  a 
free-boundary-integral  evolutionary  method  that  is  ideally  suited  for 
incompressible,  inviscid  or  nearly-inviscid  two-dimensional  flows.  The 
contours  are  the  boundaries  of  constant  density  that  are  the  sources  of 
the  flow;  e.g.,  of  constant  vorticity  regions  in  the  homogeneous  Euler 
equations  and  of  constant  mass  density  regions  in  the  stratified  Euler 
equations.  Generally,  the  velocities  of  the  contours  are  obtained  from 
integrals  (homogeneous  Euler)  or  integral  equations  (stratified  Euler, 
etc.)  on  the  contours.  Thus,  in  CD  the  evolution  of  plane  curves 
describes  the  nonlocal  and  nonlinear  dynamics  of  two-dimensional  fluid 
and  plasma  systems.  It  is  a  natural  technique  for  flows  in  unbounded 
media  because  the  Green's  function  has  a  simple  form.  After  "beyond" 
early  times, ^13,  14^  one  must  utilize  topological-change  or  "surgery" 
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algorithms  which  interconnect  and  clip  contours.  In  effect,  these 

introduce  a  smallest  scale  into  the  problem.  Dritschel  has  automated 

this  process  with  a  robust  algorithm  which  he  calls  "contour 
14 

surgery."  Although  it  violates  the  Euler  equation  conservation  laws, 
it  seems  to  give  very  good  results  "up  to"  intermediate  times  since  the 
errors  are  confined  to  the  smalls  +  scales. 

The  moment  model  is  derived  by  assuming  that  the  vortex  regions 

which  create  the  flow  are  well-separated.  A  moment  representation  for 

each  region  is  introduced  and  truncated  after  second  moments 

(elliptical  representation)  and  thus  one  obtains  two  additional  degrees 

of  freedom  for  each  centroid.  Although  an  asymptotic  description,  the 

model  works  well  for  closely  interacting  vorticies  and  resolves  the 

initial  stages  of  vortex  merger.  In  particular,  the  model  becomes 

integrable  when  applied  to  the  symmetric  merger  of  two  identical  vortex 

regions,  and  yields  explicit  necessary-and-suf f icient  conditions  for 
15 

merger.  Although  the  moment  model  is  derived  under  the  assumption  cf 
uniform  vorticity,  the  conditions  for  merger  are  in  agreement  with  our 
high-resolution  spectral  simulations  containing  initially  smooth 
circular  vortex  distributions.  This  indicates  a  certain  degree  of 
universality  of  the  merger  conditions. 

Recently  we  examined  the  problem  of  asymmetric  merger,  a  more 
realistic  problem,  and  thereby  one  richer  in  parameters.  One  may  ask: 
which  vortex  "core"  will  be  the  uictor  in  a  merger?^' ^  It  would  be 
terribly  expensive  to  answer  the  question  of  the  vortex  victor  with  a 
pseudospectral  code  and  even  with  a  contour  dynamical  code.  However, 
the  moment  model  gives  excellent  insights  into  the  victory  process  and 
yields  information  for  an  analytical  assault  on  the  unrestricted  merger 
problem. 18 
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